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We present new limits on low mass accelerator-produced 
point-like Dirac magnetic monopoles trapped and bound in 
matter surrounding the D0 collision region of the Tevatron 
at Fermilab (experiment E-882). In the context of a Drell-Yan 
mechanism, we obtain cross section limits for the production 
of monopoles with magnetic charge values of 1, 2, 3, and 6 
times the minimum Dirac charge of the order of picobarns, 
some hundred times smaller than found in similar previous 
Fermilab searches. Mass limits inferred from these cross sec- 
tion limits are presented. 

14.80.Hv,13.85.Rm,07.55.-w 



The existence of magnetic monopoles of charge g ex- 
plains the quantization of electric charge e as eg — nhc/2, 
n = ±1, ±2, ... Q|, results in the dual-symmetrization 
of Maxwell's equations 0, and is not forbidden by any 
known principles of physics. The minimum values of the 
quantization number are n = 1 according to Dirac and 
n = 2 according to Schwinger (!]] . For e being the charge 
on the electron, these values are n = 3 and n = 6, re- 
spectively, if quantization with the quark electric charges 
is allowed. Previous searches for trapped and bound 
magnetic monopoles in samples from various accelerators 
HQ, in meteorites and lunar soil |(| have been made. 
Since the Tevatron proton-antiproton (pp) collider has 
extended its integrated luminosity by a factor of some 
ten thousand over the last search of Bertani et al. [Q, 
we have taken samples exposed in the D0 experiment 
and performed a search that improves these limits. We 
use the induction method of Alvarez et al. |7j] to measure 
the magnetic charge content of macroscopic Al and Be 
samples. Indirect searches || are able to present lim- 
its beyond those energetically allowed in direct searches, 
but represent an entirely different set of assumptions as 
to their observation, and are not further discussed in 
this letter. We also do not discuss very massive GUT 
monopoles, such as those searched for in cosmic rays [||. 

This extension of limits is experimentally driven. 
Theoretical motivation, beyond the general principles 
given above, derives from the possibility that mag- 
netic monopoles generated during spontaneous symmc- 
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try breaking at the electroweak scale might give rise to 
monopoles of mass ~ 2.5 to ~ 15 TeV p0| . The search 
here can only raise the previous limits of ~ 0.1 TeV to 
~ 0.4 TeV. At the Large Hadron Collider one could ap- 
proach 2 TeV with the same techniques. 

A large warm bore cryogenic detector, similar to that of 
Longo et al. (j| , was constructed and operated at the Uni- 
versity of Oklahoma. The active elements of the detec- 
tor are two 19 cm diameter superconducting loops each 
connected to a DC SQUID (Superconducting QUantum 
Interference Device) jy]]. The Meisner effect prevents a 
change in the net flux through the loops, resulting in a 
change, or "step," in current flowing in the loops when- 
ever a magnetic charge passes through them. Samples of 
a size less than 7.5 cm in diameter by 7.5 cm in length 
are repeatedly passed through the 10 cm diameter warm 
bore centered on and perpendicular to the loops, trav- 
eling some 108 cm about the position of the supercon- 
ducting loops. In a central 64 cm region this allows for 
the magnetic effects of induced and permanent dipolc 
moments in the sample to start and return to zero on 
each up and each down traversal of the sample. High fre- 
quency noise is filtered out, while random and 1// noise 
contributions are averaged out over 20 pairs of up/down 
passes. Each up or down traversal takes some 25 sec- 
onds. A net data rate of 10 Hz is recorded for each of 
the two SQUIDs; recorded as well are the readings from 
an accelerometer, the vertical position from an optical 
encoder, the number of increments taken by the stepper 
motor, and the time. 

The SQUIDs are tuned and their transfer functions 
measured periodically according to the manufacturer's 
specifications in order to keep them operating with con- 
stant sensitivity. The absolute calibration of an expected 
signal from a Dirac monopole was made using a "pseu- 
dopole." A long thin magnetic solenoid (1.5 cm diameter 
by 100 cm length with 4710 turns per meter) carrying a 
calibrated small current gives a calculable pseudopole at 
either end. This pseudopole can either be passed through 
the warm bore of the detector in a way similar to the sam- 
ples, or it can be placed in a given position with one end 
fully extended through the SQUID loops and the solenoid 
current repeatedly switched on and off. Both methods 
were used and a value of 2.40 ± 0.04 mV/minimum Dirac 
charge was obtained. The internal random and system- 
atic local noise contributions were typically 0.2 mV, al- 
though the external systematic error of ~ 0.7 mV domi- 
nates. 
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The shape of the pseudopole curve is compared to that 
of a theoretical calculation in Fig. la and the response to 
a point magnetic dipole is compared to a narrow one from 
the data. Fig. lb shows the step from a 5.5mV (i.e., 2.3 
Dirac poles) signal from a run of a pseudopole embedded 
in an Al sample. The movement of conducting samples 
through the small residual magnetic field gradient in the 
warm bore also gives rise to an induced dipole, which 
cancels out over an up/down pair. 
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FIG. 1. "Pseudopole" curves, a) Comparison of theoretical 
monopole response to an experimental calibration and of a 
simple point dipole of one sample with that calculated from 
the theoretical response curve, b) The observed "step" for 
a pseudopole current, corresponding to 2.3 minimum Dirac 
poles, embedded in an Al sample. 

There are two kinds of samples, beryllium and alu- 
minum. A 46 cm section of the 5 cm diameter D0 
Be beam pipe, centered on the collision region, cover- 
ing nearly the full solid angle, was cut into six 7.6 cm 
pieces. Two extension cylinders of Al, each of 150 cm 
diameter by 46 cm length and of 1.26 cm thickness, also 
were cut into pieces of 7 cm by 7.6 cm or of 6 cm by 7.6 
cm, and bundled four pieces to a "sample." The effective 
solid angle acceptance Jl2| subtended by both extension 
cylinders was 0.12 of 4tt, around a cosine of the labora- 
tory polar angle of —0.82 or 0.82. There were a total of 
222 Al samples. Both the Be and Al samples were mea- 
sured using the aforementioned traversal scheme. Since a 
nylon string, which could be magnetically contaminated, 
pulled the samples through the warm bore, we also ran 
background runs between every two samples. 

The data analysis proceeded as follows. The time se- 
quences of the SQUIDs' outputs were examined inter- 
actively, and bad sections deleted pairwise (up/down); 
typically 17-18 pairs of traversals remained. A pedestal 
value (the SQUID output near the top end of each traver- 
sal) was subtracted from every voltage value along that 



up (or down) traversal. The values for each of some 90 
small ranges of vertical positions were averaged, remov- 
ing most of the random drift of the SQUIDs. The two 
SQUIDs' data were averaged, shifting one relative to the 
other by 10.1 cm in position in order to superimpose 
their dipole responses. The background samples were 
analyzed similarly and local groups of background runs 
were averaged. These backgrounds were subtracted from 
the samples' spectra. A horizontal line was fit to two 
regions, one at the lower position and one at the upper, 
as seen earlier in Fig. lb; the difference in the values of 
these two flat fits gives the step for that sample. Exam- 
ples of sample spectra are shown in Fig. 2a and 2b. The 
steps from the various samples are histogramed in Fig. 3. 
The background subtraction ensures that the distribu- 
tion of steps centers on zero, since the small effect of the 
magnetized string holding the sample has been removed. 
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FIG. 2. Sample spectra, a) Beryllium sample "SBe5P," 
and b) aluminum sample "S133A1." The observed steps are 
—0.8 mV in a) and +0.4 mV in b). The dipole signals are 
off scale in the middle regions of the plot in this vertically 
expanded view. 

The distribution of steps for the data has a mean of 
0.16 mV and an rms spread (sigma) of 0.73 mV, as shown 
in Fig. 3. One sided 90% confidence limits for monopole 
charges of n = 1 or — 1 can be obtained by considering 
the number of samples that are within 1.28 sigma of the 
n = ±1 positions, corresponding to samples outside of 
the central region of ±1.47 mV around zero. We find 8 
samples outside this central region, where 10.4 were to 
be expected from the Gaussian error. According to Feld- 
man and Cousins p3[ , the 90% confidence upper limit 
is 4.2 signal events for 8 events observed when 10 were 
expected. In order to be sure that none of the 8 outlying 
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samples were monopole candidates, we remeasured them 
along with a few control samples; all eight remeasured 
samples fell within ±1.47 mV of n = 0. 

No samples are within 1.28 sigma of the |n| > 2 po- 
sitions, the closest being 3.08 sigma away from n = —2. 
The 90% confidence upper limit is 2.4 signal events for 
zero events observed and zero expected. The upper limit 
numbers 4.2 and 2.4 are used to derive cross section up- 
per limits for \n\ — 1 and \n\ > 2, respectively. 
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FIG. 3. Histogram of steps. Vertical lines (dashed) define 
the expected positions of signals for various n. The Gaussian 
curve (dashed) corresponds to 228 measurements having an 
average value of 0.16 mV and an rms sigma of 0.73 mV. 

The acceptance of the experiment to the stopping and 
trapping of monopoles must be estimated. It is a func- 
tion of ranging out of monopoles due to the energy loss 
in bulk matter and the distribution in energy of the pro- 
duced monopoles. It is assumed that ranged out and 
stopped monopoles bind to the magnetic dipole moments 
of the appropriate nuclei, 9 Be or 27 A1 |Q. The pro- 
duction is assumed to derive from a Drell-Yan process: 
quark-antiquark annihilation to monopolc-antimonopolc 
pair via an intermediate high mass virtual photon. The 
shape of the energy distribution follows from a dimen- 
sional argument that is basic to Drell-Yan: Al 3 da/dM 
is dimensionless, where M is the invariant mass of the 
pair of monopoles. This pp cross section must include a 
threshold phase space factor and the velocity dependence 
of the monopole interaction. The threshold factor is j3, 
the velocity of the monopole in the CM system. We take 
the interaction factor to be j3 2 , since the Lorentz force for 
magnetic charges g is F = <?(H— v x D). Thus the energy 
shape of da/dM is (/3/M) 3 , convoluted with the momen- 
tum distributions of the quarks in the colliding proton 
and antiproton. In the absence of a theory of monopole 
cross sections, we are led to make these model assump- 
tions. Fig. 4a shows the shape for the cases j3 = 1 and 
(3 = v/c. The acceptance is the area of this distribution 
between two M limits divided by the area of the total 
distribution. The two M limits are M\ ow = 2m + 2Ti ow 
and Mhigh = 2m + 2Thi g h, where m is the mass of each 
produced monopole and the T's are the kinetic energies of 
monopoles which are either just entering or just ready to 



exit the sample. If the energy is too low or too high, the 
monopole will not be absorbed by that part of the detec- 
tor. The energy loss functions as summarized by Ahlen 
Jl5| were parameterized as to range, monopole mass, an- 
gle, and sample. The estimated acceptances given in Ta- 
ble I vary slowly in the mass regions of interest. 

Using the total luminosity delivered to D0, 172 ± 8 
pb _1 the number limit of monopoles, the accep- 

tance, and the solid angle coverage given above, we ob- 
tain the pp cross section limits shown in Table I. These 
limits are more than 100 times smaller than the best prior 
Tevatron limit of Bertani et al. 0] of 200 pb. 
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FIG. 4. Cross-section plots, a) da/dM for Drell-Yan pro- 
duction. The solid curve, cut between the two kinetic en- 
ergy (T) limits, yields the acceptance for stopping and trap- 
ping n — 1 monopoles in the aluminum extension cylinder. 
The dashed curve gives that for dimuons. b) The differential 
dimuon cross section is normalized to the D0 |l^] and CDF 
mji dimuon data, c) The da/dM curve of b) is multiplied by 
(r, integrated and renormalized by a factor of n 2 (137/2) 2 to 
give the cross section curves versus pair mass (2m) as shown. 
For the lunar samples, the ordinate is the proton-nucleon cross 
section. The n = 1, f3 = 1 curve is shown for comparison be- 
cause the f3 3 correction was not included in previously pub- 
lished limits. 

One can further interpret these limits as mass limits 



3 



using the scaled Drell-Yan cross sections. Here the cross 
section is taken to be n 2 (137/2) 2 larger than the Drell- 
Yan muon pair cross section, modified by /3 3 , for pp inter- 
actions measured by CDF [jl7| and by D0 . For such 
large cross sections a unitarity limit appears at an equiv- 
alent n 2 ~ 9 [[l9| . We thus use the n — 1 or 2 scalings for 
the cases n = 1 or 2, and the unitarity limit for higher 
n values in converting cross section limits into mass lim- 
its. We offer this procedure, as well as that for deter- 
mining the cross sections, as canonical. We realize that 
future theoretical work may change this interpretation 
of the data. The normalization of the Drell-Yan muon 
pair calculation (using the CTEQ5m parton distribution 
functions ^o|) to the CDF and D0 data are shown in 
Fig. 4b. The result of the integration of da/dM to a to- 
tal cross section for monopoles is shown in Fig. 4c, along 
with the upper limits for the data of Table I. The Bertani 
et al. [Q and the lunar rock sample limits of Ross et al. 
H are also shown. The corresponding interpreted lower 
mass limits are given in Table I. These limits are some 3 
times larger than those of prior searches for accelerator- 
produced monopoles trapped in matter. 
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Sample 
Afl/4n acceptance 
Mass Acceptance 
Number of Poles 
Upper limit on cross section 
Monopole Mass Limit 



Al 
0.12 
0.23 
< 4.2 
0.88 pb 
> 285 GeV 



Al 
0.12 
0.28 
< 2.4 
0.42 pb 
> 355 GeV 



Be 
0.95 
0.0065 
< 2.4 
2.3 pb 
> 325 GeV 



Be 
0.95 
0.13 
< 2.4 
0.11 pb 
> 420 GeV 



TABLE I. Acceptances, upper cross section limits, and lower mass limits, as determined in this work (at 90% CL). 
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